HIGH ENERGY NEUTRINO FLUXES FROM COSMIC ACCELERATORS 
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We constrain high-energy neutrino fluxes with the observed cosmic ray and gamma ray fluxes, 
include flavor oscillations and propagation through Earth, and show that blazars could possibly 
be detected by cubic-kilometer neutrino telescopes. 



1 Introduction 

As several cubic- kilometer neutrino telescopes are being planned, it is instructive to ask about 
the maximum event rate one should expect for these telescopes. Evidently any answer to this 
questions depends on some crucial assumptions, especially concerning the nature of possible 
neutrino sources. 

Following the discussion in Mannheim et a/.,i we derive an upper flux limit under the as- 
sumptions that (i) protons do not escape from the acceleration region (whereas neutral particles 
do), (ii) the power index of the target flux density spectrum (i.e. the synchrotron photons which 
accelerated electrons and protons can interact with photohadronically) is given by a = — 1, and 
(iii) the cosmic ray proton proton flux has the value 

Np,ohs{Ep) = 0.8 X {Ep/l GeV)-2-75 cm-^g-isr-^GeV^^ (1) 

Eq. |l| is an upper limit consistent with all observations for proton energies Ep between 3 x 10^ 
and 10^^ GeV. Assumption (ii) implies that the bolometric gamma-ray luminosity of the 
sources considered equals twice the neutrino luminosity Li,\ 
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Using particle physics, one niay obtain the muon neutrino production rate from the correspond- 
ing neutron production rate:l3 

Q^^{E) ^83.3Qn{25E). (2) 



Here and in the following particles and antiparticles are not distinguished. The cosmic ray 
proton flux (which is due to the decay of neutrons leaving the source) is given by 



Qcr{Ep) ~ Qn{Ep) X Pcsc,n{Ep), (3) 

where Pesc,n denotes the neutron escape probability. 
2 Constructing neutrino flux bounds 

The recipe for constructing a generic upper bound on the neutrino flux can now be stated as 
follows: Given some neutron production rate one may compute the corresponding production 
rates of muon neutrinos and cosmic ray protons by means of Eqs. |2| and ||. In order to compute 
the neutrino and proton fluxes / at Earth resulting from these production rates, one has to 
integrate over all sources in the universe, bearing in mind that interactions during the flight 
(e.g. Bethe-Heitler processes and pion production) change the particle spectra: 



Here z denotes the redshift factor, Vc the comoving volume, the luminosity distance, and 
dPsource/dVc the redshift distribution. M{E, z) takes into account the modification of the spec- 
trum due to interactions during the flight. For neutrinos M{E, z) = 1 holds valid. 

Finally, the expressions for I^^ and Icr thus obtained are normalized so that /cr is tangent to 
the observational cosmic ray proton flux limit (Eq. |^). If gamma-rays are overproduced for this 
normalization (i.e. if = 2Ly is greater than the bolometric diffuse gamma-ray luminosity 
the proton and neutrino fluxes are reduced appropriately. 

Fig. shows several bounds obtained by means of the formalism just described: 

1. Generic bound for optically thin sources with Qn{En) '^ffn^ exp(— £'„/£'max), the redshift 
distribution of which equals that of AGN and galaxies, u The generic bound (denoted by 
Tn-y < 1) is the envelope of the flux limits computed for each -Emax- 

2. Generic bound for optically thick sources. Here the bound (denoted by » 1) is derived 
by demanding consistency with the diffuse gamma ray background only. 

3. Generic bound for blazars. The construction parallels that of the generic bound for opti- 
cally thin sources, but we assume a spectral break between 10'' GeV and 10^^ GeV in the 
escaping cosmic ray flux due to opacity effects (still allowing for > GeV photon emission). 

4. Bounds for EGRET blazars and BL Lac objects. In computing the flux limit from Eq. |^ 
QcT and Qy are averaged over the luminosity function and redshift distribution of EGRET- 
detected blazars and BL Lacs, respectively. 

Details of the calculation of these bounds may be found in Mannheim et a/.i 



3 Event rates in neutrino telescopes 

So far, we have considered muon neutrinos only. However, as neutrinos are produced predomi- 
nantly via pion decays, one should expect the ratio Qi,^ : Qi^^ : Qur ~ 1 : 2 : for the various 
neutrino production rates. In addition, flavor oscillations during the flight to Earth lead to a 
flux ratio of ly^ : ly^ : ly^ = 1 : 1 : l.@ Hence in the following, we will assume equal numbers of 
neutrinos for the three flavors and change the upper bounds accordingly. 
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Figure 1: Left: Cosmic ray bounds on extragalactic neutrino fluxes. The generic bound for the optically thin 
{Tn-y < 1) and thick case (r„-y S> 1), the generic blazar bound (solid), the EGRET blazar bound (assuming that 
for L = 10"* erg/s the energy E-r^i at which the neutron optical depth is unity has the value i?r=i ~ 10* GeV) 
(dashed) and the BL Lac bound (assuming that Et^i ~ 10^^ GeV for L = 3 x 10** erg/s) (dot-dashed curve) are 
shown together with the Frejus limitu and the limit inferred by Waxman and Bahcall ( W&cB) with and without 
source evolution. u Right: Corresponding proton fluxes together with the cosmic ra« proton flux limit used in this 

work (straight line). Adopted from Mannheim et al.u 



Furthermore, the neutrino fluxes are altered due to neutrino-nucleon interactions during the 
crossing of the inner Earth: □ 

= -a.Mm.iE) + Y.f dE'^iE', E)U,(E') + decay of r. 

k ^ 

In this equation, i and k denote neutrino flavors, t the column number density. dak^i{E\ E) jdE 
constitutes the differential cross section for turning a neutrino of flavor k and energy E' into one 
of flavor i and energy E. (Xj^tot is the total (charged and neutral current) cross section of fj. 

Limiting ourselves to muon neutrinos, we may write the differential muon event rate dN /dE 
in a neutrino telescope with an effective area ^eff as a sum of two terms covering the muons due 
to neutrino-nucleon interactions inside and outside the detector: 

Here dacc/dE constitutes the differential cross section for charged current interactions, p the 
density of the detector medium (i.e. p = 1 g/cm^ for water), L the detector length, and rup the 
proton mass. Eq{E,x) denotes the energy a muon must have so that after crossing a distance x 
of the detector medium (and thus suffering radiative losses) it retains the energy E. 

The event rates corresponding to various neutrino flux bounds are shown in Fig. ^ for 
horizontal and vertical incidence, assuming a cubic kilometer telescope. 

4 Discussion 



Fig. H clearly shows that whereas EGRET-blazar-like neutrino spectra might imply an event 
rate of about ten events per year and steradian, generic blazar models could allow neutrino 
event rates of several hundred events per year and steradian (most sources would then not be 




Figure 2: Left: Muon event rates in a cubic kilometer Cerenkov telescope, assuming flavor oscillation and hori- 
zontal incidence (no propagation through the Earth). Energy losses outside the detector are taken into account. 
The curves correspond to the respective curves of Fig. 0. Right: Corresponding rates for vertical incidence 
(propagation along the diameter of the Earth). The generic blazar event rate lies below 1 year~^sr~^. 

gamma ray emitters above GeV). This exceeds the hmit obtained by Waxman and Bahcall @ 
for optically thin sources (which is comparable to our BL Lac case) by one to two orders of 
magnitude, reflecting the different assumptions. 

Any derivation of an upper neutrino flux limit faces two important unknowns: □ Firstly, 
intergalactic magnetic fields (which affect protons but have no influence on neutrinos) might 
considerably change the number of neutrinos compatible with cosmic ray observations. Secondly, 
the same applies to exotic processes such as the decay of superheavy particles. 

The planned cubic-kilometer neutrino telescopes will reach the sensitivity necessary to probe 
the assumptions that enter the flux-limit calculations. If photohadronic processes play an im- 
portant role in blazars, we have shown that they will very likely be detected by such telescopes. 
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